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ABSTRACT 

We report on the first wide-field, high-precision astrometric analysis of the 13 extragalactic radio sources of the complete S5 polar cap 
sample at 15.4 GHz. We describe new algorithms developed to enable the use of differenced phase delays in wide-field astrometric 
observations and discuss the impact of using differenced phase delays on the precision of the wide-field astrometric analysis. From this 
global fit, we obtained estimates of the relative source positions with precisions ranging from 14 to 200 /jas at 15.4 GHz, depending 
on the angular separation of the sources (from ~1.6 to ~20. 8 degrees). These precisions are ~10 times higher than the achievable 
precisions using the phase-reference mapping technique. 

Key words. Astrometry - techniques:interferometric - galaxies:quasars:general - galaxies:BL Lacertae objects:general - radio con- 
tinuum: general 



1. Introduction 

In the past few years, we have carried out a series of very-long- 
baseline-interferometry (VLBI) observations, using the Very 
Long Baseline Array (VLBA) at 8.4, 15.4, and 43 GHz, aimed at 
studying the absolute kinematics of a complete sample of extra- 
galactic radio sources using astrometric techniques. The target 
of our programme is the "complete S5 polar cap sample", con- 
sisting of 13 radio sources from the S5 survey (Kiihr et al. 1981 , 
Eckart et al. 119861) . All sources in this sample have flux densities 
higher than 0.2 Jy at the epoch of our observations and well- 
defined International Celestial Reference Frame (ICRF-Ext.2) 
positions (see Fey et al. 2004*). The relative source separations, 
less than about 15° for neighbouring sources, should allow for 
typical astrometric precisions in the range of 20 to 100 //as, de- 
pending on the observing frequency. (Lower frequencies trans- 
late into lower precisions.) Even though the essence of our global 
differenced phase delay astrometry is the same as that of phase- 
reference astrometry (e. g., Beasley & Conway 1995), there are 
important differences between them. 

On the one hand, the coordinates of one source (the target 
source) are determined in phase-reference astrometry with re- 
spect to the coordinates of another source (the reference source). 
In our global analysis, we use data from all the 13 sources of 
the S5 polar cap sample simultaneously, in a unique fit, which 
increases the precision of the astrometry by a factor ~10, as we 
will see in Sect. 4. On the other hand, phase-reference astrom- 
etry requires small angular separations between the sources (up 
to a few degrees). In our global analysis, this requirement is dra- 
matically relaxed. Consequently, we should be able to apply a 
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global differenced phase delay astrometric analysis, similar to 
the analysis reported here, to a set of sources spread across the 
whole sky. 

We have already described the goals of our astrometric pro- 
gramme in Ros et al. (2001 hereafter Paper I) and presented 
VLBA maps obtained at two different epochs at 8.4 GHz in 
Paper I and at 15.4 GHz in Perez-Torres et al. (2004, hereafter 
Paper II). In this paper, we present the first conclusive phase 
delay astrometric results of this programme. We describe in 
Sect. 3 the details of our astrometric analysis and our algo- 
rithm developed in-house from which we automatically derive 
the 2n- ambiguities inherent in the phase delay observable. We 
also discuss the contribution of differenced observables in the as- 
trometric precision for sources distributed across a large portion 
of the sky. 



2. Observations and maps 

We observed the complete S5 polar cap sample at 15.4 GHz on 
15 June 2000 using all the antennas of the VLBA for 24 hr. The 
observations took place in subsets of 3 or 4 radio sources, which 
were always observed at high elevations. The sources of each 
subset were observed cyclically for about 2 hr. On-source scans 
were 60 sec long, with a small time gap (10-20 sec) to slew the 
VLBA antennas. Thus, one complete observing cycle was about 
5 min long. 

This observing mode resulted in a total observation time 
for each radio source of about 4hr (see Fig. []]). Data were 
cross-correlated at the Array Operation Center of the National 
Radio Astronomy Observatory (NRAO). We used the NRAO 
Astronomical Image Processing System (aips) for the calibration 
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Fig. 1. Sky distribution of the polar cap sample sources and 
schematic representation of the schedule of our observations. 
Each circle points to the sources included in the same observ- 
ing block, typically lasting 2-4 hours with a duty cycle of about 
5 minutes. The numbers associated to the circles represent the 
chronological order of the observation blocks. 



of the visibilities. We aligned the visibility phases through the 
whole frequency band (for all sources and times) by first fringe- 
fitting the single-band delays of one scan of source 1803+784 
and then applying the estimated corrections to all the visibilities. 
Thus, another fringe-fitting using the multi-band delays provided 
the new phase corrections for all the observations. The visibility 
amplitude calibration was performed using the system tempera- 
tures and gain curves from each antenna. For imaging, we trans- 
ferred the data into the program difmap (Shepherd et al. 1 19751 1 
and made several iterations of phase and gain self-calibration 
until obtaining high-quality images (with residuals close to the 
thermal noise). The images of all the sources corresponding to 
this epoch were analysed and published in Paper II. The struc- 
tures of the S5 polar cap sources for our epoch, and other epochs 
at 15.4 GHz and 8.4 GHz, were discussed in Papers II and I, re- 
spectively. 

3. Astrometric analysis 

We performed a differenced astrometry analysis based on the 
phase-connection of the phase delays (e.g. Shapiro et al. 119791 
Guirado et al. 119951 Ros et al. [1999 and Perez-Torres et al. 
2000), but with some substantial differences. For the astrometric 
fits (i.e., estimates of the coordinates of the sources via weighted 
least-square fits of the astrometric observables), we used the 
University of Valencia Precision Astrometry Package (uvpap), an 
extensively improved version of the vlbi3 program (Robertson, 
119751 see also Ros et al. 1999 for details of the vlbi3 model). The 
main improvements of uvpap permit the use of the JPL ephemeris 
binary tables and compute the relativistic effects of the Solar 
System bodies using the Consensus Model (McCarthy & Petit 
2003). We outline the main steps in the astrometric analysis. 



We used aips to obtain the group delay, phase delay, and rate 
from each observation of each of the thirteen sources, af- 
ter subtracting all the contributions from the structure of the 
sources, thus referring the phase delays to the phase centre of 
the images. We discarded the data from the antennas Mauna 
Kea and St. Croix, because of the high noise in their corre- 
sponding fringes. 

2. We predicted the number of phase cycles between consecu- 
tive observations of each of the thirteen sources to permit us 
to "connect" the phase delays (e. g., Shapiro et al. 1 1 979l >; the 
computation of the number of phase cycles was performed 
by comparing the phase delays with the modelled delays ob- 
tained from uvpap, using a fit of the clock drifts of the VLBA 
antennas and the tropospheric zenith delays to the group de- 
lay and rate data. 

3. We refined this phase delay connection using an algo- 
rithm that imposes the nullity of all the closure phases (see 
Sect. 3.1 for details). 

4. We computed the differenced phase delays among the 
sources observed in the same blocks (see Fig. [T}. These dif- 
ferenced delays are largely free of unmodelled effects from 
the troposphere, ionosphere, and antenna electronics (e. g., 
Marcaide et al. |1994| l. 

5. We estimated the positions of the 13 sources of the S5 polar 
cap sample via a global weighted least-square analysis of the 
undifferenced and differenced data. 

In the last step, we fitted the positions of all the sources with 
respect to the phase centre of 0454+844 (whose coordinates 
were fixed parameters in our fit) along with the tropospheric 
zenith delays (see below) and clock drifts for each antenna. We 
chose the source 0454+844 as the reference source, not only be- 
cause of its position (roughly in the geometrical centre of the sky 
distribution of the sample, minimising the sum of distances to the 
other sources), but also because it was the most observed source 
in the schedule, being part of most of the observation blocks (see 
Fig. [J}. Taking this source as reference therefore provides more 
stability to the astrometric fit. 

Regarding the propagation medium, we modelled the tropo- 
spheric zenith delay at each station as a piecewise-linear func- 
tion, with a separation of 6 hours between nodes, that is, there 
were five nodes for each antenna. A priori values at the nodes 
were calculated from local surface temperature, pressure, and 
humidity measured at each station, based on the Saastamoinen 
model (Saastamoinen I1973K We used the dry Chao mapping 
function (Chao et al. |1974| l to determine the tropospheric de- 
lays at other elevations than the zenith. Changing the mapping 
function did not alter our results much; for instance, the astro- 
metric corrections of the relative source positions (see below) 
obtained with the use of the Global Mapping Function (Boehm 
et al. 120061 1 differed less than 8 //as with respect to those ob- 
tained using the Chao mapping function. For the ionosphere, 
we used global ionospheric maps at the epoch of our observa- 
tions derived from GPS data and generated by the Center for 
Orbit Determination in Europe (CODE). These maps (in IONEX 
format) provide the Earth's total electron content (TEC), which 
we transformed into delays over each station using the atps task 
tecor (e. g., Walker & Chatteriee [T999l Ros et al. |2"0"00j i. 

3.1. Phase-connection algorithm 

We used the astrometric model in uvpap to predict the num- 
ber of phase cycles (27r-ambiguities) between consecutive ob- 
servations of the same source by means of the residual delay 
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Fig. 2. Distribution of the residual delay rates for all the base- 
lines, sources, and scans of our observations. 



rate. Since the average time separation between observations is 
~ 180 sec and the phase cycle at 15.4 GHz is ~65ps, the residual 
rates have to be lower than (65ps/180sec)=0.36ps/sec to en- 
sure a good phase connection. In Fig. [2] we plot the distribution 
of the residual delay rates of our observations. Residual delay 
rates with absolute values higher than 0.36ps/sec in the distri- 
bution, most of them corresponding to the longest baselines and 
the weakest sources, are not uncommon. Thus, this preliminary 
phase connection is far from perfect and both time- and baseline- 
dependent phase cycles remain uncorrected in our data. In previ- 
ous works, with less antennas and less sources (e. g., Guirado et 
al. 119951 Ros et al. II 9991 ), these additional cycles could be cor- 
rected by inspection of the phase residuals. However, this pro- 
cedure is unmanageable with the amount of data in our obser- 
vations. Therefore, we have devised an algorithm to automate 
the correction of these unmodelled phase cycles. The algorithm 
works as follows: 

1 . For each source and observing time, we compute the closure 
phases corresponding to all triplets of antennas. From the 
non-zero closure phases, we determine the baseline involved 
most frequently. 

2. We perform an ambiguity check on this baseline. By ambigu- 
ity check, we mean shifts of the phase delays of that baseline 
by adding or subtracting one phase delay cycle. We perform 
a positive (i. e., addition of a cycle) and a negative (i. e., sub- 
traction of a cycle) shift. We then compute the scores corre- 
sponding to both shifts. By score, we mean the number of 
closure phases approaching zero minus the number of clo- 
sure phases distancing from zero after the shift. 

3. We select the shift with highest score and modify the data 
with such a shift. 

4. We iterate this process (1 to 3). We keep checking on the 
ambiguities, until all the closure phases (i. e., closure phase 
delays) are made zero for the source and time selected. 

We apply this algorithm to all scans in our observations. 
Actually, to ease the work, for consecutive scans of the same 
source, we apply the corrections found in the previous scan be- 
fore performing the ambiguity check. We have tested this "auto- 
matic connector" with synthetic data for different scenarios and 
have obtained excellent results. An example scenario, with a re- 
markably high noise level, consists of delays equal to a random 
number of cycles (up to 5) added to randomly selected baselines 



every random number of scans (with an average of 50 scans be- 
tween random cycles) in a dataset of 1000 scans and 10 antennas. 
Under such circumstances, the automatic connector finds all the 
random baseline-dependent cycles without introducing changes 
in the antenna-based overall constant cycles (see below). We re- 
peated this test several times (with a different number of anten- 
nas and even in worse noise conditions) and the connector never 
introduced changes in the antenna-based overall cycles. 

3.2. Antenna-based ambiguities 

The antenna-based ambiguities are offsets of the phase delay, 
consisting of a given integer number of phase delay cycles that 
depend on each antenna and source. These ambiguities do not 
affect the phase closures and, thus, are completely transparent 
to the automatic connection algorithm described above. These 
antenna-based cycles appear very clearly in the residuals of the 
differenced delay observables, but can also be detected in the 
undifferenced residuals. To correct these antenna-based ambigu- 
ities, we applied another algorithm based on a smoothness cri- 
terion, which analyses variations between differenced (and also 
undifferenced) residuals of neighbouring scans that are, in mod- 
ulus, close to, or larger than, a phase cycle. For each scan, the 
algorithm 

1 . finds these variations. 

2. analyses whether these jumps have an antenna-based struc- 
ture. 

3. corrects the antenna-based phase cycle in the observations. 

This algorithm was applied in a "bootstrapping" manner, 
i. e., beginning with a subset of 3 close-by antennas (Kitt Peak, 
Pie Town, and Los Alamos) and adding more antennas (one at 
a time) when all the residuals of the subset of antennas were 
finally smoothed. We show an example of one iteration of this 
smoothness criterion algorithm in Fig. [3] 

3.3. Overall (source-based) ambiguities 

Since we had phase-connected the data for each source indepen- 
dently, we still had to determine, for each antenna, the overall 
source-based ambiguity, that is, the integer number of phase cy- 
cles by which the phase delay of one source is offset from the 
others. This offset cannot be totally absorbed by atmospheric, 
clock, or astrometric corrections, and it can notably affect the 
astrometric results at our precision level. 

We determined the overall ambiguities by again following 
an iterative process: first, we estimated the overall ambiguities 
in our weighted least-square fit. The overall ambiguities clos- 
est to an integer number of cycles of phase delay were set to be 
exactly equal to that integer number of cycles of phase delay. 
Then, we repeated the astrometric fit to obtain new estimates of 
the remaining overall ambiguities. Progressively, all the ambi- 
guities were fixed to integer numbers of cycles of phase delay. 
For the complete set of overall ambiguities, the maximum devi- 
ation with respect to an integer number of cycles turned out to 
be less than one fourth of a cycle of phase delay (in fact, 40% 
of all the overall ambiguities deviated less than one tenth from 
their closest phase cycle integers). 

3.4. Differenced observables in the global fit 

The use of differenced observables corresponding to a given pair 
of sources will increase the precision in the determination of the 
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Fig. 3. Example of the effect of an antenna-dependent unmodelled phase cycle in the differenced observables associated to some 
baselines and source pairs (see source codes in Tabled. The baseline and source pair corresponding to each figure are indicated in 
the upper-left corner (BR refers to Brewster, OV to Owens Valley, and NL to North Liberty). The dashed lines represent the delays 
+0. 06515ns, corresponding to +1 and -1 cycle of phase at the observing frequency. The smoothness criterion algorithm (see Sect. 
3.2) infers, from (a) and (b), that source 04 is the one affected by the unmodelled phase cycle and, from (c), that the affected antenna 
is OV; (d) corroborates the deductions extracted from the other figures. 



relative positions of such pair (i. e., the position of one of the 
sources with respect to the position of the other). Differenced 
observables for a total of 24 source pairs (see Table [3]l could 
be computed. This "network" of differenced observations intro- 
duces redundancies for sources that appear in more than one 
pair. These constraints in the degrees of freedom for the po- 
sitions allow for an increase in the precision, not only of the 
relative source positions, but also of their absolute coordinates, 
although to a lower degree. The advantage of the redundancy 
introduced by the network of differenced observations can only 
be used when more than one pair of sources is available. Thus, 
the use of all our data in a unique fit provides more robust re- 
sults than the sub-division of the observations in individual sets 
of source pairs fitted separately. 



Unlike other differenced analyses, in which one source of the 
pairs was always fixed in the fit, in our global scheme we have 
differenced observables constructed with pairs of sources whose 
coordinates are being simultaneously estimated in the astromet- 
ric fit (except the coordinates of source 04, which are kept fixed 
in the fit). This led us to reconsider (see Appendix A) the concept 



of "changes in the relative position" of a pair of sources (a, b), 
which is now defined as 



Aa r a = Aat, - Aa a - A6 a sin(a/, - a a ) tan((5/,) 
A5 r a el = A5 h - A6 a cos(a b - a a ) 



where (Aa r a , AS r a ) is the change in the relative position of 
source b with respect to source a, (a a ,6 a ) and (at, 6b) are the 
ICRF-Ext.2 right ascensions and declinations of sources a and 
b, respectively, and {Aa a , A6 a ) and (Aa/,, A6h) are the astromet- 
ric corrections for sources a and b resulting from our fit. These 
expressions have been obtained considering the curvature of the 
celestial sphere when none of the two sources is kept fixed in 
our fit. See Appendix A for more details. We notice, in partic- 
ular, that Aa"' and A8 r f l need not be the negative of Aa r f and 
A6 r a el , respectively, since the different orientations of the sources 
in the sky will produce a combination of their respective changes 
in right ascension and declination. This is particularly true for 
the sources close to the celestial pole, as it is the case for the 
high-declination sources of the S5 polar cap sample. 
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Table 1. Fixed parameters in our astrometric fit and their uncer- 
tainties, according to IERS. 



200 



Fixed parameter 



A-priori uncertainty 



0454+844 position 
Earth Pole 
Site Coordinates 
UT1 - UTC 



300 yuas (in a and S) 
0.7 mas (in a and 5) 
2 cm in each coordinate (x,y,z) 
0.04 ms 



3.5. Global phase delay astrometric fit 

We use both differenced and undifferenced observations in the 
same (global) fit, and the latter are included for the fit to remain 
sensitive to antenna-dependent parameters (i. e., clock drifts and 
zenith delays). We scale the standard deviations of the differ- 
enced and undifferenced phase delays separately in such a way 
that, for each type of data and for each baseline and source, the 
rms of the postfit residuals is unity. The ratio between standard 
deviations of the differenced and undifferenced observations is 
-0.67. 

To estimate the astrometric uncertainties, we allowed for 
variations in the site coordinates, the coordinates of the Earth's 
pole, and UT1-UTC in an auxiliary fit, but with their adjust- 
ments constrained by their a priori values and their standard de- 
viations (see Table [TJ through the use of an a priori covariance 
matrix. We then used the final covariance matrix of this auxil- 
iary fit to estimate the final astrometric uncertainties, which now 
include all the contributions and correlations between the param- 
eters of the geometry of the interferometer and the propagation 
medium. We also scaled the uncertainties of all the fitted param- 
eters to make the reduced^ 2 equal to unity. 

Regarding the tropospheric zenith delays, the standard devi- 
ations of the fitted nodes of the piece-wise linear functions of our 
model were ~0.01 ns. These are smaller than the expected uncer- 
tainties due to random variations in the wet component of the tro- 
pospheric delay at each site (~0.1 ns between nodes; Treuhaft & 
Lanyi 1987). These a priori, extra random variations are not con- 
templated in the computation of our astrometric uncertainties 
because such unmodelled variations constitute a large fraction of 
the final rms (see figure|U), which indirectly affects the astromet- 
ric uncertainties after their scaling to obtain a reduced y 1 equal 
to unity. 

Our global astrometric fit provides the standard deviations of 
the absolute positions of the sources. Given that we use Eq.Q]to 
calculate the changes in the relative positions of the sources, the 
corresponding standard deviations for these changes are given 
by the following expression: 



Ay, 



-^-C(x- x) 



[ cr(A6 r a el ) = 



(2) 



where jc,- is the set of corrections (Actb, A6b, Aa a , A6 a ), and 
C(Xj, Xj) is the element (i, j) of the covariance matrix. 

We tested the robustness of our astrometric results by repeat- 
ing the fit with different reference sources. The new absolute 
source positions were compatible with those from the first fit at 
1 cr level. The relative source positions were compatible at 2 cr 
level. 



1 Fixing the uncertainties of the nodes of the tropospheric delay mod- 
els to 0. 1 ns, would result in an increase in the astrometric uncertainties 
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Fig. 5. Uncertainties in the relative coordinates a (triangles) and 
6 (squares) of all the pairs of sources as a function of their sepa- 
rations. 



4. Results and discussion 

We present the results of our astrometric analysis in Tables|2]and 
|3] Table [2] shows the astrometric corrections and corresponding 
uncertainties to the absolute source positions given by ICRF- 
Ext.2. Table [3] shows: i) the changes in the relative coordinates 
of the 24 source pairs, computed using Eq. [TJ ii) the angular 
separations of the source pairs, computed from their ICRF-Ext.2 
coordinates; and iii) our estimated corrections to the source pair 
angular separations. 

As an example of the quality of our fit, we show in Fig. 
|4]the residuals of the undifferenced and differenced phase de- 
lays corresponding to all the observed sources and one of the 
longest baselines (Hancock - Kitt Peak). Note the cancella- 
tion of systematic effects in the differenced data, which are still 
noticeable in the undifferenced data (probably unmodelled at- 
mospheric effects). The rms of the undifferenced delays for all 
sources and baselines range between 55 ps (baseline Fort Davis 
- North Liberty observing the source 02) and 6 ps (baseline Kitt 
Peak - Pie Town observing the source 00). The rms of the dif- 
ferenced delays range from 35 ps (baseline Hancock - Owens 
Valley observing the pair 11-04) to 2.2 ps (baseline Brewster - 
Los Alamos observing the pair 18-17). The latter rms is some- 
what smaller. 

We emphasise that the relative coordinates of each pair in 
Table [3] are not the simple subtractions of the absolute coordi- 
nates of the sources forming such a pair in Table [2] (except for 
those pairs with reference to the source 04; see Eq.QJi. Also, the 
uncertainties in the relative coordinates are much smaller than 
those of the absolute coordinates. This is a consequence of, first, 
the natural cancellation of systematic errors in the differenced 
observables and, second, the correlations between the absolute 
positions estimated in the fit that account for all posible global 
shifts of the sources of each pair. 

As could be expected, the standard deviations of the absolute 
positions are roughly the same for all sources; however, since we 
are using differenced observables in our fit, the standard devia- 
tions of the relative positions of the 24 pairs are strongly depen- 
dent on the separation of the sources that form the pairs. This 
effect can be clearly seen in Figs.[5]and[6] Figure|5]shows the un- 
certainties in the relative a and S for the 24 pairs as a function of 
their separations. This behaviour is roughly linear, as predicted 
by the empirical formulae given by Shapiro et al. (119791 ) and 
corroborated by simulations of astrometric VLBA observations 
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performed by Pradel, Chariot & Lestrade (2006). This linearity 
is more clearly seen in Fig. [6] which shows the uncertainties of 
the estimates of the separations as a function of such separations. 
However, our astrometric uncertainties are ~10 times smaller 
than those predicted by Shapiro et al. ( 1 19791 ) and those esti- 
mated by Pradel, Chariot & Lestrade (2006) from their simula- 
tions, for typical declinations of the S5 polar cap sample sources. 
This improvement in precision is probably due to the fact that 
we have analysed the 13 sources simultaneously, instead of us- 
ing the usual two-source (target/reference) scheme. Therefore, 
our global analysis brings more constraints on the position of 
the sources, which produces an immediate benefit in the preci- 
sion of the astrometry for particularly high-declination sources. 
Moreover, unlike the phase-reference mapping technique, the 
tropospheric delay model is re-estimated along with the source 



positions in our least-square fit, which also contributes to min- 
imise the^- 2 . 

4. 1 . Brief discussion of some selected sources 

All the astrometric corrections shown in Tables [2] and [3] are 
corrections to the ICRF-Ext.2 positions. These ICRF-Ext.2 po- 
sitions, determined using the group delay observable obtained 
from dual frequency observations at 8.4 & 2.3 GHz, are not di- 
rectly comparable to our 15.4 GHz phase delay position esti- 
mates: first, the ICRF-Ext.2 positions are not defined with re- 
spect to any specific phase centre on the source maps, which 
our positions are; second, even if the source structure is not as- 
trometrically significant, source opacity effects could be present 
while comparing the source positions determined at 8.4, 2.3, and 
15.4 GHz. 

From Table [2] we can calculate mean corrections to the 
ICRF-Ext.2 coordinates. We find mean corrections of 278 and 
170 juas in right ascension and declination, respectively. Our val- 
ues are similar to those found by Jacobs et al. (2004). These 
authors are pursuing the extension of the ICRF-Ext.2 to 24 (K- 
band) and 43 GHz (Q-band). From VLBA global observations 
they found an agreement between the X/S-frame and K-band 
frame to within 330 and 590 //as in right ascension and decli- 
nation, respectively, comparable to our mean corrections. 

Structure effects must contribute to the corrections of the 
separations of the source pairs given in Table [3] Only seven out 
of the thirteen S5 sources are defining ICRF-Ext.2 sources and 
only one of them (source 07) is unresolved (following the defini- 
tion used by Chariot [19901 ) at X-band. Actually, a detailed look 
at the structure of the radio sources can provide some hints to ex- 
plain the relatively large corrections in the last column of Table 
[3] For instance, pairs with source 01 as member, which has an 
X-band structure index of 4 (very extended, see Chariot 1990) 
show large corrections to the ICRF-Ext.2 coordinates; likewise, 
the corrections corresponding to pairs containing source 1 1 can 
be explained in terms of a southwest bending of the jet near the 
core (see maps in Paper I & II). However, a detailed interpreta- 
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Table 2. List of the S5 polar cap sample sources. 



Source 


Alias 


Flux density 6 


ocr_ 


ICRF-Ext.2 position (J2000.0) 


Astrometric corrections'' 






(mJvi _ 


(mJy/beam) 


(Y 




Air ( usis.} 


AS (ua^ 


0016+731 


00 


823.1 


0.5 


nnh iqITI jc 7Q(S491 S 


73° 27' 30 01750" 


_C20 + 320 


-560 + 200 


01 53+744 


01 


339.7 


0.7 


01 n 57 m 34 964908 s 


74° 42' 43 22998" 


-1 50 + 300 


-810 + 180 

O ± \J _!_ luU 


0212+735 


02 


2519.1 


1.2 


nnh 17m in (n 3373s 


73° 49' 32 62179" 


_f,20 + 310 


-280 + 190 


0454+844 


04 


225.3 


0.6 


05 n 08 m 42 363503 s 


84° 32' 04 54402" 


± 300 


± 300 


061 5+820 


06 


445.7 


0.7 


06 n 26 m 03 0061 88 s 

\J\J £J\J \JJ .\J\)\J 1 OO 


82° 02' 25 56764" 


110+ 1 60 


200 + 1 00 

L.\)\J i 1 \J\J 


0716+714 


07 


1063.4 


0.8 


07 h 53.448459 s 


71° 20' 36.36339" 


220 ± 360 


-100 ± 210 


0836+710 


08 


1772.9 


0.8 


08 h 41 m 24.365262 s 


70° 53' 42.17301" 


280 ± 360 


-400 ± 220 


1039+811 


10 


848.4 


1.0 


10 h 44 m 23.062554 s 


80° 54' 39.44303" 


660 + 190 


-30 ± 110 


1150+812 


11 


988.3 


1.0 


ll h 53 m 12.499130 s 


80° 58' 29.15451" 


-450 + 180 


-640 + 120 


1749+701 


17 


468.7 


0.8 


17 h 48 m 32.840231 s 


70° 05' 50.76882" 


460 ± 380 


140 ± 230 


1803+784 


18 


2334.5 


0.7 


18 h 00 m 45.683914 s 


78° 28' 04.01849" 


70 ± 230 


130 ± 140 


1928+738 


19 


3088.6 


1.1 


19 h 27 m 48.495167 s 


73° 58' 01.57010" 


90 + 310 


0+ 190 


2007+777 


20 


1271.0 


0.5 


20 h 05 m 30.998511 s 


77° 52' 43.24763" 


400 ± 240 


-250 + 150 



" Aliases of source names used throughout the paper. 

* Source flux densities at the epoch of our observations, obtained from hybrid mapping using natural weighting. 

3 times the root-mean-square of the residual images. By residual image, we mean the Fourier transform of the difference between measured 
visibilities and model visibilities. 

d Astrometric corrections to the ICRF-Ext.2 source positions, resulting from our work. Source 04, fixed at the ICRF-Ext.2 position, is taken as 
reference. 



Table 3. Astrometric results for all the observed source pairs. 



Source pair^ 


Change in the relative coordinates^ 


ICRF-Ext.2 angular separation 


Correction to angular separation 




Aa rel (//as) 


AS rel (pas) 


(deg) 


yuas 


01-00 


-990 + 70 


180 ±40 


6.770731378 


850 ± 70 


01-02 


-390 + 16 


520 ± 14 


1.614936659 


-620 ± 20 


04-01 


-150 + 110 


-810 ± 140 


12.284604495 


810 ± 140 


04-02 


-620 + 110 


-280 + 140 


12.703755162 


450 ± 140 


04-06 


110 + 40 


200 + 40 


3.332785200 


-110 ±40 


04-07 


220 ± 150 


-100 ± 160 


14.393224831 


150 ± 170 


04-10 


660 ± 100 


-30+ 110 


10.085113231 


390 ± 120 


06-07 


-80 + 90 


-300 + 120 


11.093608166 


280 ± 130 


06-10 


350 ± 90 


-120 + 80 


9.134095175 


350 ± 100 


08-04 


-400 ± 180 


240 ± 80 


16.400909237 


460 ± 190 


08-07 


-180 + 70 


270 ± 20 


6.419140221 


150 ± 70 


08-10 


730 ± 110 


310 ±90 


12.179115583 


760 ± 140 


11-04 


-350 ± 120 


-120 ± 100 


11.407147545 


350 ± 140 


11-10 


930 ± 40 


580 ± 20 


2.699376639 


-1020 ±30 


11-18 


1280 ± 130 


110 ± 140 


14.839207812 


700 ± 170 


11-20 


1530 + 130 


-600 ± 200 


18.626089045 


1180 ± 210 


18-00 


-1050 + 160 


-550 ± 200 


20.844902015 


-140 ±240 


18-04 


-60 + 80 


130 ± 190 


16.901504141 


-140 ± 190 


18-17 


340 + 50 


10 ± 100 


8.408293010 


-30 ± 100 


18-20 


260 + 70 


-360 + 40 


6.342305345 


370 + 70 


19-17 


350 + 70 


140 + 70 


8.491093690 


-360 ± 90 


19-20 


330 ± 30 


-250 ± 50 


4.521891420 


-15 ± 50 


20-00 


-1150+ 130 


-450 ± 140 


15.450912042 


-500 ± 180 


20-02 


-910 + 150 


-300 + 200 


20.618479511 


-300 + 240 



" See aliases defined in Table|2] 

b Defined as the change in coordinates of the second source with respect to the first one, according to Eq. IA.5l The uncertainties are estimated 
using Eq.|2] 



tion of the astrometric information obtained by comparing with 
the results from other epochs and frequencies in our astrometric 
project is beyond the scope of this paper and will be published 
elsewhere. 



5. Conclusions 

We report here on the first high-precision, wide-field astromet- 
ric results at 15.4 GHz of our multi-frequency monitoring of 
the S5 polar cap sample. To obtain those results we first de- 
veloped the package uvpap, an extensively improved version of 
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the vlbi3 program (Robertson 119751 ). The ability of uvpap to use 
differenced phase delays, along with newly developed phase- 
connection algorithms, enables us the use of differenced phase 
delays in global astrometric observations. We discuss the im- 
pact of the differenced phase delays on our global astrometric 
analysis and show that their use increases the precision of the 
relative source positions by a factor of ~10 compared with the 
precision achievable using the phase-reference technique with a 
pair of sources (Pradel, Chariot & Lestrade 2006). The astromet- 
ric precisions obtained linearly decrease (from 14 to 200 //as) 
as the separations between the sources increase (from ~1.6 to 
~20.8 degrees), with the result that the fractional errors in deter- 
mining the separations of all the studied source pairs are similar 
(~3xl0- 9 ). 

We obtain some large corrections for the relative coordinates 
and separations of the sources. From all the 24 pairs studied, 10 
have separation corrections above 500 /ias and, of those, 4 have 
corrections above 900 /ias. These corrections could be caused by 
opacity effects (our observations are at 15.4 GHz, and the ICRF- 
Ext.2 positions are based on 8.4 and 2.3 GHz observations) and 
by source structure effects (we are relating our astrometric posi- 
tions to the phase centres of the source maps, that is, the peaks 
of brightness; the ICRF-Ext.2 positions are not well-defined on 
the source structures). Other wide-field, high-precision astro- 
metric analyses of these sources at other frequencies and other 
epochs are currently under way. This multi-epoch and multi- 
frequency study will eventually provide spectral information and 
the absolute kinematics for all sources in the sample. Ultimately, 
we expect to provide a definitive test of the stationarity of the 
innermost radio-source cores, associated to the massive black 
holes, which is a basic tenet of the standard jet interaction model 
(Blandford & Konigl [T979l ). In addition, our results will be an 
excellent complement to future yuas-precise astrometry at optical 
wavelengths. 
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Appendix A: Relative position changes for 
free-moving sources in the sky 

Let a a and 5 a be the right ascension and declination of source 
a, and at, and 5b the right ascension and declination of source b. 
Then, the relative position of source b with respect to a is 



a h 
8 b ■ 



■ a a 
5„. 



(A.l) 



In the simple case that one of the two sources (the refer- 
ence source a) is kept fixed in the sky, the change in the relative 
position between this pair of sources can be well defined from 
equations (A.l): 



Aa r a el = Aab 
Adf = A5 b 



(A.2) 



where a is the reference (fixed) source and b the target (free) 
source; Aab and A5h are the corrections to the position of source 
b, maintaining the source a fixed in the fit. However, Eqs. (A.l) 
are not appropriate when the two sources can change their po- 
sitions in the fit, since the curvature of the celestial coordinate 
system affects the robustness of Eqs. (A.l) under a global shift 
of the source pair. To illustrate with a simple example for two 
sources separated by 12 hours in right ascension, a global shift, 
e, in declination would originate a change of 2e in the relative 
declination of these sources, according to (A.l). (The declina- 
tions of these sources will change with opposite signs.) The rela- 
tive coordinates between this pair of sources would, then, appear 
to change dramatically under a global shift of the pair. In other 
words, we need to define a reference point in the sky to measure 
the shift of b with respect to a. In our analysis, for each pair of 
sources, we select the nominal position of the reference source 
as the reference point for the study of that particular pair. In prac- 
tice, such a selection is equivalent to applying a global rotation 
% in such a way that the source a is rotated back to its initial 
(i. e. nominal) position. According to this rotation, the change in 
the coordinates of source a is 



a a + Aa a -> a a 
5 a + A5 a 5 a 



(A3) 



where Aa a and A5 a are the corrections to the right ascension and 
declination of source a. This rotation, ^R, will have the corre- 
sponding effect on the coordinates of source b, which are 



a h + Aa b 
5 b + A5 b 



a h + Aal' b + Aa, 



5 b + A5>;; b + A5 b 



(AA) 



where Aa" a ' b and A5" a ' b are the changes in right ascension and dec- 
lination that the rotation % causes on the position of source b. 
Thus, the change in the relative coordinates of b with respect to 
a will be 



Aa r f = Aa b + Aa'" 



A5 r f = A5 h + A5 



(A.5) 
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Fig. A.l. Graphic representation of the rotation % that we apply to sources a and b. In (a.l), the coordinates of both sources, a and 
b, have been corrected an amount Aa = (Aa a , A6 a ) and Ab = {Aa/,, Ad/,), respectively. In (a.2), we apply a rotation % that brings 
source a back to its ICRF position (maintaining P ang constant) and causes the shift A m ft = (Ao^, AS'" b ) on source b. The total shift 
of b will therefore be the addition of the corrections Ab and A" 



ab 



When the coordinates of both sources, a and b, are corrected 
in the astrometric fit, we must use these equations instead of Eq. 

From all possible rotations % in the sky that move the source 
a back to its nominal position, we selected the one that causes 
the direction of the arc between a and b (i. e., the position angle 
of b with respect to a) to remain unchanged. According to the 
sine theorem (see Fig. |A. Il l: 

cos 5b sin (a a - ab) = sin A sin P ang (A. 6) 

where A is the arclength between a and b, and P ang is the position 
angle of b with respect to a. In our case, this can be expressed as 



cos (5i, + A5b) sin (arj - a a + Aa b - Aa a ) = 
cos (6 b + A5" a ' b + A5 h ) sin (a b - a a + Aa™ h + Aa b ), (A.l) 

For a first-order approximation of the astrometric corrections 
(~0.5 mas, see Tabled, the condition lA.7l together with the con- 
stancy of the arclength between a and b under the rotation 9v, is 
satisfied if, and only if 

Aa™ b = -Aa a - A5 a sinfe - a a ) tan(^) (A.8) 
A<5™ = -A6 a cos(a h - a a ). 

When we apply these relationships to Eq. IA.51 we obtain Eq. 
Q]of this paper directly: 

Aa r a el = Aa b - Aa a - A6 a sin(a£ - a a ) tan(6b) (A.9) 
A5"' = A5 b - A5 a cos(a b - a a ). 



